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We calculate, in a model, the beam spin asymmetry in semi-inclusive jet production in deep 
inelastic scattering. This twist-3, T-odd observable is non-zero due to final state strong interactions. 
With reasonable choices for the parameters, one finds an asymmetry of several percent, about the 
size seen experimentally. We present the result both as an explicit asymmetry calculation and as a 
model calculation of the new transverse-momentum dependent distribution function g^. 
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I. INTRODUCTION 

Semi-inclusive deep inelastic scattering (SIDIS) pro- 
vides a way to reach a more detailed understanding of 
the structure of a hadronic target [U, 0, [1] . The informa- 
tion gained can be codified into transverse momentum 
dependent distribution functions (TMD's), which in lead- 
ing twist are related to quark probability distributions. 
Single spin asymmetries provide a method to isolate and 
measure specific TMD's often those of higher twist. 

Experimentally, both beam and target spin asymmetries 
have been observed 0] ■ 

This paper is concerned with calculations of beam sin- 
gle spin asymmetries in semi-inclusive production of a 
quark jet. A beam spin asymmetry can be viewed in a 
number of ways. At the simplest level, finding a beam 
spin asymmetry means defining "up" and "down," say by 
taking a cross product of the incoming electron spin 
with the virtual photon three-momentum q, and then see- 
ing if the outgoing hadron or jet is more likely to emerge 
with transverse momentum up or down relative to the 
electron scattering plane as just defined. Going further, 
a jet is a collimated spray of hadrons, where all the final- 
hadrons in the jet are summed and their phase space 
integrated, and we represent the jet momentum by the 
quark momentum pq. A more formal observable that 

corresponds to the beam spin asymmetry is Se x q ■ Pq- 
Non-zero beam spin asymmetry means a non-zero aver- 
age, weighted by the cross section, value of this observ- 
able. One can define an azimuthal angle (f> as the angle 
between the electron scattering plane and the plane de- 
fined by the photon and quark momenta. The observ- 
able has a sine/) dependence, so that a non-zero beam 
spin asymmetry requires a term in the differential cross 
section proportional to sin (f>, and the relative size of this 
term is the measure of the beam spin asymmetry that we 
shall use in the body of this paper. ^ 

The electron interacts with the hadrons by photon cx- 
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^ If parity is violated, as in neutrino scattering, there is an addi- 
tional term in the single-spin asymmetry proportional to sin 2(/) 



change, so one may ask how the photon obtains the spin 
information and carries it across. The detail is that lon- 
gitudinally polarized electrons, which we consider here, 
give rise to photons that have a vector polarization S^, 
and that this vector polarization is not parallel to the 
photon momentum although it is in the electron scat- 
tering plane. (Virtual photons emitted by unpolarized 
electrons have only tensor polarization.) An alternative 
beam spin asymmetry observable is x q- pg. 

The observable Se x q ■ pq is odd under time-reversal 
(T-odd). Hence the beam spin asymmetry is and must be 
zero in a lowest order calculation. However, it is not zero 
in higher orders, as hadronic final state interactions lead 
to a relative phase between the longitudinal and trans- 
verse amplitudes, and it is interference between these two 
amplitudes gives the beam spin asymmetry. 

Single spin asymmetries in semi-inclusive reactions 
have their own history. Part of the history is common 
to beam and target single spin asymmetries, since both 
are T-odd, though it developed first in the target spin 
asymmetry context. Early on Sivers Q pointed out that 
target spin asymmetry could be used to measure a new 
distribution function for the proton. However, it was 
quickly noticed that the asymmetry was T-odd and so it 
was first thought that the asymmetry had to be zero [911 . 
So it remained until Brodsky, Hwang, and Schmidt [lO| 
demonstrated explicitly in a model the statement already 
made above, that because of final state interactions the 
T-odd operators could give experimentally non-zero ex- 
pectation values. Shortly thereafter, this statement was 
confirmed from other viewpoints [ill, ■ The same sim- 
ple model has since been applied to a number of other 
target spin dependent TMD's 

The other part of the history impacted upon questions 
of factorization in SIDIS. It is believed or hoped — and 
perhaps proofs are becoming available [3, HE] for lead- 
ing twist — that cross sections for SIDIS can be written 
as a sum of terms, where each term is a convolution 
of a distribution function that depends only on the tar- 
get and a fragmentation function that depends only on 
the final state. For jets observed in the final state, the 
fragmentation function would be either absent or triv- 
ial. However, until relatively recently, the lists of possi- 
ble target TMD's contained none that could give a beam 
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spin asymmetry when the observed hadron was the whole 
quark jet 0, A problem arose when we, in an earlier 
note and Metz and Schlegel directly calculated 
the beam spin asymmetry and found a non-zero result 
without using any T-odd fragmentation function. One 
could then speculate that factorization did not work for 
this situation since there was then no parton distribution 
associated with this asymmetry. However, Bacchetta et 
al. [l^ and Goeke et al. [13] , have made progress by un- 
covering an initial state distribution function g-'-, which 
was originally dismissed because it was T-odd, and which 
precisely gives a beam spin asymmetry with an outgoing 
jet. Factorization for observable matching g-^ ^ and twist- 
3 distributions, has yet to be proved. (We might note 
that if specific hadrons rather than a jet are observed, 
there are and were known mechanisms to obtain a beam 
spin asymmetry [1, [^, [2ll, Also, there were early 

perturbative QCD calculations including gluon loop dia- 
grams which showed beam asymmetries of the order of a 
percent [23|.) 



The model we use [10|, [17|, |18|, l2J] is a simple one where 
the proton is represented as a bound state of a quark and 
a scalar diquark, and where the hadronic final state in- 
teraction is produced by a gluon exchange. The beam 
spin asymmetry involves both longitudinal and trans- 
verse photon polarizations. The longitudinal matrix el- 
ement (for spin-1/2 quarks) is subleading for large Q. 
For reasons of electromagnetic gauge invariance, we keep 
all Feynman diagrams where the photon attaches to a 
charged particle. Hence we have three diagrams, in gen- 
eral, in lowest order, and there are also three diagrams 
in one- loop order that have an imaginary part. The dia- 
grams are shown below in Figs. [Hand [21 For calculations 
where only the transverse amplitudes are needed to lead- 
ing order in l/Q, such as the unpolarized cross section or 
the target spin asymmetry [3, only the diagrams where 
the photon interacts directly with the quark are needed. 



We show our calculations in Sec. [TTl with the results 
for the beam spin asymmetry given both as explicit for- 
mulas, for the limit that the masses and jet transverse 
momentum are all smaller than the invariant momentum 
transfer Q, and as plots for selected values of the parame- 
ters. Our calculation of the beam single spin asymmetry 
does not depend on factorization theorems. Now, with 
modern knowledge, finding a beam spin asymmetry does 
not disprove factorization, and if factorization is correct 
for the situation at hand, wc can interpret the beam spin 
asymmetry in terms of the distribution function . We 
do so in Sec. IHIl and include a few further comments 
about factorization. The beam spin asymmetry we dis- 
cuss here has no mention of hadronic spin, so one expects 
we could exhibit a similar asymmetry in a model with 
scalar quarks. We remark further about this in Sec. IIVI 
Finally, we end with a summary in Sec. FVl 



II. CALCULATIONS 

We present calculations of the beam single spin asym- 
metry in a model where the nucleon is represented as a 
bound state of a spin-1/2 quark q and a scalar diquark 5, 
and where the momentum of the quark in the final state 
is measured. Physically, the latter means that the quark 
jet in the final state is observed, and one measures its to- 
tal momentum, but not any detailed features that would 
give information about, for example, the polarization of 
the quark. The fundamentals of the model are the same 
as in Ref. (Toj . 

The process, with momenta and helicities shown, is 

e{l, K) + N{p, A) -> e{l', h,) + g(pi, A') + S{p2) . (1) 

We treat the electron as massless, so that electron helic- 
ity hf. = ±1/2 is conserved. The masses of the nucleon, 
quark, and scalar diquark will be M, m, and tos, respec- 
tively. The general matrix element is 

M = ^(e{l',h,)\jMl,he)) 

X (q{pi,X')S{p2)\r\N{p,X)), (2) 

where ji, and J'^ represent electromagnetic current oper- 
ators for electrons and hadrons, respectively, and the last 
matrix element will often be abbreviated J'^(A, A'). 

The five-fold differential cross section for a polarized 
electron but unpolarized hadron is 



da 



(3) 

where Ee and E'^ are the electron energies in the target 
rest frame ("lab"), q = I - V , ^ -g^ = [p + qf , 
and ^ is the quark three-momentum magnitude in the 
hadronic center-of-mass. The lepton tensor is 

L^^ihe) = (e{l,he)\jf_t\e{l\he))(e{l',he)\j„\e{l,he)) 

= Ll, + {2K)L^,, (4) 

where "5"' or '!A" indicate terms that are symmetric or 
antisymmetric in the indices {fJ.,v). The hadron tensor 
is here just 



(5) 



A, A' 



Experimenters measure a single spin asymmetry de- 
fined from 



da cc ( 1 + (2/ie)^i(7 ' 



whence 



A^'°'^sin< 



(6) 



(7) 
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The pointed brackets indicate an averaging over the 
quark azimuthal angle (p. (Such an averaging was also 
done for the unpolarized term in Eq. 

After some modest effort, the last equation can be con- 
verted to 



L2 



V2e(l -e)lmW 

' eWLL + \{Wii+W22) 



(8) 



The indices "L" on the hadronic tensor indicate contrac- 
tion with the exchanged photon's longitudinal polariza- 
tion vector e^, 



L2 



W 

(9) 

Plain e is the photon polarization parameter, given by 

(10) 
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where v ^ E(. — E'^ and 9^ is the electron scattering angle 
in the target rest frame. An alternative expression for e 
using 



IS 



l-y-yV(4r) 
l-2/ + 2/V2 + 2/V(4r) 



(11) 



(12) 



In the scaling limit (y ^ co with y fixed), r — *■ oo. 

For spin-1/2 quarks, the transverse amplitudes domi- 
nate the longitudinal ones by a power of Q in the scaling 
limit, and e goes to a constant in the same limit. Hence 
the transverse terms arc dominant in the denominator of 
Eq. ([8]) above, and the asymmetry J^^lf" falls like l/Q at 
high Q. 

We need to make some comments about gauge invari- 
ance. Such discussion was unnecessary in [lO[ . Those au- 
thors calculated a target spin asymmetry which depends 
only on amplitudes with transversely polarized photons, 
which automatically conserve the electromagnetic cur- 
rent. Beam spin asymmetry depends on an interference 
between longitudinal and transverse photon amplitudes, 
and longitudinal amplitudes are gauge sensitive. 

Gauge invariance requires that we allow the photon to 
interact with all charged particles in the Feynman dia- 
gram. For 7* + ^ g + S", this gives in lowest order 
the diagrams in Fig. [TJ There is a "quark graph," a "di- 
quark graph," and a "proton pole graph." In [lO|, only 
the quark graph [(a,0)] was included. In the large Q 
limit and for transverse photons, the other two graphs 
are smaller by at least a factor l/Q. For longitudinal 
photons, the additional graphs are not suppressed. 

In general, one has all three diagrams. Charge conser- 
vation. Cat = ei +62, ensures gauge invariance. (Charges 
Cat, ei, and 62, arc for the nuclcon, quark, and diquark, 
respectively.) 




(a,0) 




(b,0) 




(c,0) 

FIG. 1: Lowest order graphs. 



We calculate in a Breit frame, and find it convenient 
to adopt much of the notation of Ref. [l^- We use the 
notation for the outgoing quark transverse momen- 
tum. The external momenta, in light front coordinates 
((j± = ± g3), are 

9 = 9", = (-Q,Q,0_l) , 



P = 



Pi = 



X ' Q 



P2 - Q 



l^x A l+m^ X Al+ml 



The electron scattering plane defines the x-z plane, with 
the outgoing electron having a positive component in the 
X direction. Momenta p and q are exact. In the high 
limit X becomes the Bjorken variable [xBj = Q^/'^P • <?]; 
momenta pi, 2 are accurate to 0{1/Q). 

In the Breit frame and using gauge invariance, one can 
show 

{eL), J^ = J+. (14) 

The zero-loop diagrams give an amplitude J^'^^(A, A') 
which is the sum of 



u(pi,A')7^(l*i- fi+m)u{p,\), 



(pi - qf - m2 
•^r^^o) " ~ {P2 - 7?- ml ^^^^ ~ '^^^ ^(Pi^^'MP' ^) : 
^(%) = - w^l\p ^ AO (^+ ^ + M) r<P, A) . 



(15) 



Thus 
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where, 



•^(a,0)('^''^') 
'^(t,0)('^'^') 



a;(l 



eig 1- X Q_ 

l2 + A2 x/^ 



m2 + 



2eig 1 — a; 
2-a; 1 



623- 



2 1 



/X 



TO 

X 



1-X 



iX5 



A, A' 



(m + xil/)|A_L|e-''^'^ - TO|A_L|e'' 



X4> 



2 + Ai <5A,A'(™ + a:A/)-A,5A,-A'|A±| 



m^A.A' +A,5A,-A'|A±|e*^'^ 



(16) 
(17) 
(18) 

(19) 



r 



(20) 



Matrix elements J^^^ and J^^ are subleading in 
l/Q- Regarding >/(aO)' only need to know that 
k(U(^'^')l = I^K0)(A,A')|. 

The calculated beam spin asymmetry, Eq. ([8]), is zero 
if we have only the lowest order amplitudes. To obtain a 
non-zero beam spin asymmetry we need to include next- 



to-lcading order amplitudes and obtain a non-zero phase 
relative to the lowest order terms. The diagrams are 
shown in Fig [21 Only one-loop diagrams that give an 
imaginary part are shown. 

In the model, the final state interaction is mediated by 
a gluon that couples to a strong charge which is carried by 
the quark (strong charge Cg) and diquark (strong charge 
(— e^)) but not by the proton. One obtains the QCD 
equivalent by letting Cf (47rQ;s) with Cp = 4/3. 

The diagrams give 



-ieigel 



d^k u{pi,y) i^i + 2?^2- 1^) (^ + m) 7^ (^- ^ + to) u{p, A) 



167r4 [k^ - m2 -1 


-^e] [{k-qr- 


rn? -f 


- le] [[p + 


q — ky — rrig -{- ie] [{k — piY - 


+ ie] ' 


d'^k 


u{p^,x'){A + 


2^2 


-^){}t + 


to) u{p, A) {2p^' + q^' - 2k^') 




167r4 [k^ - m2 -t 


- if\ [{p — kY — 


m|H 


' le] [{p + 


q — fc)2 — ?T!,g + ie] [{k — pi)^ - 


- + ze] 


d^k 


u{pi,\') {A 


+ 2j 




+ to)(^+^+A/)7^u(p, A) 






-f ie] [{p + qf - 


hp 


+ ie] [{P ^ 


~ q — fc)2 — m| + ie] [(/c — pi)^ 


- + ie] 



jU ■ 2 

"^(c,i) — -''■'^Ngsg I 7734-77;^ ™2 , ,-,1 r/„ , „\2 » r2 , r/„ , „ rvi ™2 , rT2 ..2 , ' y'^^) 

where fi in the denominator is a small gluon mass temporarily included. While present in the amplitudes, it will not 
appear in the final answer for beam spin asymmetry Ai^jj indicating that the considered observable is infrared-safe. 
The absorptive part (the part that is imaginary relative to the lowest order) is obtained by substituting 

k^-L^+^eip + q-kr-ml + .e = '^"^''^^'^ " "^'^^^"^ ^ ^ ' ' '^'^^ ' ^''^ 

As a result, the 4-dimensional integration over the loop momentum in Eq. (|2ip is reduced to a 2-dimensional angular 
integration, making the result also safe from ultraviolet divergence. In the numerical calculations, we found projections 
of the electromagnetic current Eq. pT|) on six independent helicity amplitudes of the process j* + N ^ q + S, and 
Alu was expressed in terms of these helicity amplitudes. Two-dimensional integration with respect to polar and 
azimuthal angles of the intermediate quark was done numerically with Mathematica in the 7* -I- A'^ center-of-mass 
frame. Independence of Alu on the gluon cut-off mass /z was verified, and providing a cross-check of the numerical 
calculation. The plots of the asymmetry Am shown later in this paper were obtained using the procedure described 
in this paragraph. 

We may also obtain approximate analytic formulas for the beam spin asymmetries in the limit where the masses 
and transverse jet momentum are much smaller than the invariant momentum transfer Q, by keeping leading terms 
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in a 1/Q expansion. To begin, wc have 



AbsJ(';,)(A,A') 
AbsJ(';^)(A,A') 
AbsJ(';,^(A,A') 



87r2 Q4 

e2gel 1 
87r2 l-2;g4 



Z(pi, A') ^2 (/^ + m) A) (2p ^ + g'^ - 2fc^) 
(fc_i-Ai)V/i^ 



BNgel a;2 1 /■^2^ u(pi, A') ^2 (/^ + m) (:^+ ^ + M) 7''w(p, A) 



87r2 (1-.t)2Q4 



(23) 



The fcx integrals have a finite upper limit because the quark and diquark in the intermediate state are now on-shell. 
Doing the integrals yields 



AbsJ(2,,i)(A,A') = ^jf,,o)(A,A')iA^ 



eigei l- x 



8^(m2 + Ai) ,/^2 + a2 



X |iA(5A,A 
AbsJ(+i)(A,A') = |^J+o)(A,A')i^A, 



~ 2 

i\<t> _|_ „i\<t> 



(m + xiV/)|Aj.|e-''-^*' + ^— e' 
|A±' 



>LAm, (24) 



2ei5e2 1 — x 1 



87r(TO2 _^ A2_ 



22 + A2 



+ A5a^_a 



(m + a-M) ((m^ + A^) - (m^ + A^)) La™ + ^(m^ + A^) (Lqa - 1) 



-(m^ + Ai) + f 1 - (m^ + Ai)^ La™ - (m^ + Ai)(LQA - 1) 



|Ax|e 



iXtp 



AbsJ+,)(A,A') = ||j+„)(A,A')|iA,. + iQA 
Abs J+i)(A,A') = ||j+o)(A,A')|lap + i:QA - l| 



2(1 -:g) 
2 - a: f ' 



(25) 

(26) 
(27) 



or 



Abs J+ = Abs J+ + Abs J+ + Abs J+ = 



1-x 1 



'(a,l) 



(M) 



(Cl) 



X <^ ,5 



A, A' 



87r(m2 + A^) 

2ei{m + xM) {m^ - m^) L Am - ( (2ei + 62) Xqa - 2 (d + 63) )m(m2 + A^) 



m2 + A2_ 



A5a.- 



-w? - xmM — ^— ^ La™ + (2ei + 62) Lqa - 2 (ei + 62) (?«^ + A^) 

IAi|2/ 1 — xV / 



(28) 



using charge conservation. 
Here 



Laiu = In- 



Lau = In ■ 



Lqa = In 



QHl-x) 
{fh? + 



The general result for the beam spin asymmetry is 



Ai |Aj 



47r 



Ai 



(772 + A/a;)2 + A2_ Q 

?fi'^ 1 In 



(29) 



(30) 



,r2 2 2 2ei -a;e2 .2 
MX — m + ™ 



2ei(l - a;) 



7772 + A2 2ei + 62 a; , Q'^il-x) 61+62 x 
' -in 



2ei 1 — a; 



(?7i2 + K\ 



61 \ — X 



3o{k-p, 




< 



pq 
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0.14; 
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0.08 ; 
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0.04; 
0.02; 

0; 
-0.02 ; 



• CLAS 2003 

up-quarkjet 

down-quark jet 

quark+diquark pole 

— quark+proton pole 
quark-only 



= 4.25 GeV 

beam 




0.05 0.1 0.15 0.2 



0.25 
X 



0.3 0.35 0.4 0.45 0.5 



Bj 



FIG. 3: Beam single spin asymmetry for quark jets of different 
flavors and beam energy of 4.25 GeV. Data are from Ref. 01. 



FIG. 2: One-loop graphs that have imaginary parts. 



The asymmetry falls like l/Q, as expected for a twist-3 
observable. 

For numerical work we need to make choices for the 
charges, and we will show results for several possibilities. 
Two choices, will be for up-quark and down-quark jets 
from a proton, where ei is the charge for the stated quark, 
ejv = 1; and 62 is the remainder. Another proton- like 
choice has just the nucleon and quark diagrams, with 
Gat = ei and with the diquark electrically neutral. The 
explicit asymmetry for this case is. 



J sin (f 



A 




1 



In- 



I 



1 -a; 



In- 



A2 



1 



(31) 



A further case keeps just the quark and diquark di- 
agrams. This case is neutron-like, since the nucleon is 
chosen electrically neutral, and 62 = — ei. This is the 
prescription used by Metz and Schlegel [l^ to ensure 
gauge invariance, and it gives 



A 



LU 



4tt 



rri^ -I- A^ 



lA, 




{m + Mxf + Q 



■In 



2(1- 
QHl~x) 



■ fri' 



In- 



(m2 -h A2_) 



(32) 



The result can be checked against [T^ in the limit to ^ 0. 

Finally, in our earlier note [13] , we implemented gauge 
invariance using only the quark diagram [(a)], but mod- 
ifing its amplitude by subtracting out the scalar 



PQ 



0.16; 
0.14; 
0.12; 
0.1; 
0.08 ; 
0.06; 
0.04; 
0.02; 

0; 
-0.02 ; 



CLAS 
up-quarkjet 
down-quark jet 
quark-Fdiquark pole 
quark-Fproton pole 
quark-only 



= 5.7 GeV 

beam 




^ , \ X \ , \ X , \ X \ , \ X I , d 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
^Bj 



FIG. 4: Beam single spin asymmetry for quark jets of different 
flavors and beam energy of 5.7 GeV. Data are from Ref. Q]. 
Notation is as in Fig[3] 



part [2£ 



^ -q^ '^ 

r 



(33) 



For comparison, we report the result from this "g^ sub- 
traction" method here. It is 



iv? + A^ 



47r 



(m + Ma;)2 + Ai Q 




l-2x 



In 



2(1-2;) 
Q2(i_a;) l-2a; 



2(1 - x) ^^2 + Ai) X 2(1 - x) 



(34) 



Numerical results are given in the form of plots of A^^^ 
vs. Xbj for three different beam energies. For all the plots. 
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we use parameter values M = 939 MeV, m = 300 MeV, 
mg = 800 MeV, and as = 0.3. The quark transverse 
momentum |Aj^| is fixed at 0.4 GeV. Results for beam 
energies of 4.25 GeV, 5.7 GeV, and 27.5 GeV, for quark 
jets of different flavors and with different versions of the 
model, are shown in Figs.[3l|4l and[5l The dashed-double- 
dotted curve (up to a corrected sign) matches results of 
our earlier report p7| . For lower beam energies, we start 
the plots from larger values of XBj to ensure the condi- 
tion > holds. The solid and dashed curves, re- 
spectively, describe the asymmetry in electroproduction 
of up- and down-quark jets if all three diagrams of Fig. 
2 are included. If the diquark is assumed to be not in- 
teracting with electromagnetic probe, only the diagrams 
Fig. 2(a) and (c) contribute, and the corresponding curve 
is dashed-dotted. If the nucleon target is electrically neu- 
tral, then the virtual photon couples only to the quark 
and an (oppositely-charged) diquark; the asymmetry ap- 
pears equal for both quark jet flavors. This situation is 
described by a dotted curve in Figs. [31 IH andO 

The results on these plots are from the numerically 
integrated results obtained as described just after Eq. 
([22)) . The analytic forms valid for large Q are essentially 
identical to these results for the Ee = 27.5 GeV case, and 
are relatively 10-20% lower for the other two energies. 

Remembering that our calculations refer to quark jets 
with no fragmentation into hadrons, we can see that both 
the sign and the magnitude of the asymmetry Alu com- 
pare favorably with experimental data from HERMES 
and JLab CLAS. We also note that the experiment ob- 
served somewhat stronger supression of the asymmetry 
at larger xbj compared to our model predictions. 



0.1 



D 

< 0.08 
S 0.06 
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in 
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■ 7t , HERMES prelim. 

up-quark jet 

do wn-quark jet 

L quark+diquark pole 
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E^, =27.5 GeV - 
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' 'i' 1 


" 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



FIG. 5: Beam single spin asymmetry for quark jets of dif- 
ferent flavors and beam energy of 27.5 GeV. Also shown are 
preliminary data from HERMES [^. Notation is as in FigO 



dent (TMD) distribution function g-^{x, A^), for the sit- 
uation that only the outgoing quark jet is observed and 
in the limit of neglecting the quark mass. Until relatively 
recently, for the quark jet case, the known lists of initial 
state quark distribution functions contained none that 
could give a beam spin asymmetry, although beam spin 
asymmetry was known in cases that involved quark frag- 
mentation functions that gave information about quark 
spin correlations. Now, with the work of YuanQ, Bac- 
chetta et al. [l^, and Goeke et al. the distribution 
functions hi{x,p\) and g^{x,p'^) have been uncovered, 
and they allow an observable beam spin asymmetry with- 
out any real or implied spin measurement in the final 
state. The asymmetry due to is proportional to the 
quark mass, and we consider in this Section only the 
massless quark case, and for simplicity also let ms = M. 

From Eqs. (15-17) of [3, suitably tuned for the 
present case and using the sign convention from ^20^ (that 
is opposite to [l^), one obtains 



V2e(l - e) 



\A^\xg^{x,A] 



fi{x 



Al) 



(35) 



The TMD function /i(a;,Al) gives the probability dis- 
tribution of quarks in a hadron, and is related to the 
hadron matrix by 



2{2Trf2Miyfi{x,A] 



(36) 



Thus, in the present model. 



fiix,Al) = 



1 g'a-xr 

167r3 {m^ + Aiy 



{(TO + a;A/)2 + Ai|, 
(37) 



(before m is set to zero). We obtain a value for g from 
the normalization condition. 



iV, = 1 = / dx 



d^A^hix,Al). 



(38) 



The distribution fi is calculated from lowest order di- 
agrams, so that the momentum A^ is identical to the 
transverse momentum of the quark in the initial state. 
We pick a cutoff pertinent to the process at hand, in 
particular we choose it equal to the maximum value of 
transverse momentum for a real quark in the final state. 



A2 = A2 



= Q'- 



1-x 
ix 



(39) 



Thus, 



III. THE DISTRIBUTION FUNCTION 

The general expression for the beam spin asymmetry 
can be summarized by the transverse momentum dcpen- 



167r2 



31 

18 



(40) 
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The result for 5"*" is 




The quantity in curly brackets is the same as Eq. (j30p . 
except that the mass m is set to zero both there and in 
Eq. (|^D|) . If desired, one may also substitute the con- 
tents of the curly brackets for the special cases shown in 
Eqs. dSI]), dSH), or (111. 

Again, our calculation of the beam single spin asym- 
metry does not depend on factorization. Obtaining 
does depend on the analyses of Refs. and [l^l, and 
hence does assume factorization. Factorization appears 
to have been proved for some twist-2 transverse momen- 
tum dependent distributions with transversely polarized 
targets [13, EE]- For longitudinal polarization and for 
twist-3 distributions, the question of factorization has 
been raised far in the past Q and continues to be dis- 
cussed {e.g., H^l). Factorization for these cases remains 
assumed rather than proved, and this should be borne in 
mind. 

There is a log term in g^ , mirroring a similar term 
in the result for the beam spin asymmetry. The possi- 
bility of such terms arises because the asymmetry and 
g^ are non-zero only when loop corrections are included. 
Loops, however, do not always lead to logQ^ terms; it 
depends on the details of the integrals. For contrast, 
there is the twist-2 target spin asymmetry or Sivers effect, 
where the integrals are insensitive to their upper limits in 
a similar one-loop calculation [To| and consequently one 
finds no logQ^ term. We interpret the logQ^ simply as 
5"*" possessing a logarithmic scaling violation within this 
model. It is not, as far as we can see, related to questions 
regarding factorization. 

A plot of ^-'-(a;, A^) is shown in Fig. [5] for |A_l| = 0.4 
GeV and varying x, for two of the special cases described 
earlier. Except for the logarithmic dependence in- 
duced in the normalization condition, the result is inde- 
pendent of beam energy. (The range of x for which the 
curve is plotted is fixed on the left by requiring to 
be above AP and on the right by requiring that W , the 
final state hadronic mass, is above 2 GeV. This figure as 
shown is for Ebeam = 27.5 GeV.) 



IV. BEAM SPIN ASYMMETRY FOR SCALAR 
QUARKS 

The fact that there is a beam single-spin asymmetry 
for an outgoing jet does not depend on the spin of the 
target or of its constituents. We emphasize that here 
by calculating for hypothetical scalar quark bound 



1.2 



^ 0.8 



0.6 - 



"cn 0.4 



0.2 - quark-i-diquark pole 

quark-i-proton pole 

I ' ' ' ' 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
X 

FIG. 6: The distribution function g^{x,A]_) for |Ax| = 0.4 
GeV. The two special cases are described in detail in the text. 
A short summary is that the nucleon in the quark-|-diqurk pole 
case is overall electrically neutral, and in the quark-fproton 
pole case has unit charge. In both cases all flavors of quark 
in the final state are summed. 



states. 

Some of the details are different, however. In particu- 
lar, for scalars the longitudinal amplitude dominates the 
transverse one at high Q, in contrast to the spinor case. 
Also, considering just the case that the diquark is neu- 
tral, only the quark pole diagrams are needed for the 
cross section and asymmetry calculation. To wit, using 
notation from Fig. [1] we have 



eig 



[pl ^ 9)^ ~ 
{p + qY - AP 



(2pi - qf , 
{2p + qr. (42) 



Both graphs, with ejv = ei, are needed to satisfy gauge 
invariance. However, for the longitudinal amplitudes. 



'(a,0) 



J7 



ei.g(2 - x) 



(43) 



(^■0) (l-x)Q ' 

so that the proton pole graph contribution is subleading 
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at high Q. Similarly for the transverse amplitudes, 

eig{l - x) 



J, 



(a,0)J_ 



(c,0)± 



Oi 



(44) 



After also calculating the relevant one-loop graphs, the 
beam single-spin asymmetry for semi-inclusive jet pro- 
duction is 



_ 4a, ^2e{l-e)Q + 

^L(7 ~ ~ ' ~ ~ 



A2 



2Al 



|Ax| 



(45) 

for Q large compared to masses and to |A^|. For the 
curious, evaluating the above scalar formula numerically 
with the same values of masses and |A^| mentioned pre- 
viously gives about the level of the (presumably spin-1/2 
quark) data in Fig. [3] if we use a rather small a = 0.075. 



V. DISCUSSION 

We have calculated, in a definite model, the beam spin 
asymmetry for semi-inclusive inelastic electron scatter- 
ing. Experimentally, beam spin asymmetry has been 
seen at JLab/CLAS and at Hermes. Our particular cal- 
culation considered beam spin dependent asymmetry in 
jet production, where all the final-hadrons in the jet are 
summed and their phase space integrated. 

Beam spin asymmetry means an asymmetry in the di- 
rection of an outgoing hadron relative to a direction set 
by the polarization of the incoming beam. In the case 
at hand, define an observable by first defining a normal 
from the electron spin and the photon momentum, Se x q. 
Representing the momentum of the jet by the momentum 
Pq of the quark, the beam spin asymmetry observable is 
Se X q ■ Pq- If one defines a coordinate system, this ob- 
servable is proportional to sine/), where (j) is the azimuthal 
angle between the lepton scattering and jet production 
planes. Having a beam spin asymmetry means that the 
expectation value, weighted by the cross section, of this 
observable is non-zero. 

As a useful detail, the electron spin direction is com- 
municated to the hadrons through the photon, and if the 
electrons arc polarized, the photon aquires a vector po- 
larization S-y which is in the electron scattering plane but 



not parallel to q. An alternative beam spin asymmetry 
observable is x q ■ pq. The observable is odd under 
time reversal. Further, this observable is reminiscent of 
the single spin asymmetry observable for a polarized pro- 
ton target, x q ■ Pq, which measures the Sivers func- 
tion [10, l2g and which also was thought for some time 
to be absent or zero because of time reversal invariance. 

The beam spin asymmetry is zero in a lowest order 
calculation, as expected for a T-odd observable. It is 
not zero in higher orders. Hadronic final state inter- 
actions lead to a relative phase between the longitudi- 
nal and transverse amplitudes, and interference between 
these two amplitudes gives the beam spin asymmetry. 

The model in which we calculated the beam spin asym- 
metry is a simple one where the proton is represented 
as a bound state of a quark and a scalar diquark, and 
where the hadronic final state interaction is produced by 
a gluon exchange. Gauge invariance requires that we in- 
clude all diagrams where the photon couples to a charged 
hadron. There are three diagrams in lowest order, and 
also three relevant diagrams at one-loop order. Results 
for the beam spin asymmetry A''^'^' have been given in 
the body of the paper both numerically in plots and as 
formulas for the case when masses and transverse jet mo- 
mentum are much smaller than the invariant momentum 
transfer Q. 

Quark spin not crucial to obtaining a beam spin asym- 
metry. The qualitative effect persists in an equivalent 
model with scalar quarks, although the numerical results 
are not the same. 

In modern times, the expression for the beam spin 
asymmetry in jet production can in general be given in 
terms of the initial state transverse momentum depen- 
dent parton distribution called g'^{x,k'^), where k± is 
the transverse momentum of the struck quark relative 
to the parent proton. The model thus gives an explicit 
result for , which is shown in Eq. (j4ip . 
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